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■Abstract 



Ch "In this paper we discuss the normal and superconducting state properties of two pnictide superconductors, LaOFeAs and LaONiAs, 
^ wsing Migdal-EIiashberg theory and density functional perturbation theory. For pure LaOFeAs, the calculated electron-phonon 
coupUng constant A = 0.21 and logarithmic-averaged frequency = 206^, give a maximum Tc of 0.8 K, using the standard 



Migdal-Eliashberg theory. Inclusion of multiband effects increases the Tc only marginally. To reproduce the experimental Tc, a 
^"^"5-6 times larger coupling constant would be needed. Our results indicate that standard electron-phonon coupling is not sufficient to 
C/3 explain superconductivity in the whole family of Fe-As based superconductors. At the same time, the electron-phonon coupling in 
Ni-As based compounds is much stronger and its normal and superconducting state properties can be well described by standard 
Migdal-Eliashberg theory. 
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I Introduction 

Strong-coupling electron-phonon (EP) theory, also known as 
Migdal-Eliashberg (ME) theory, was developed in the 60 's and 
70's to describe the physical properties of superconducting ele- 
mental metals and alloys, which could not be described by the 
weak-coupling BCS approach. 

, The electronic and phononic systems are described by a set 
t)f coupled diagrammatic equations (ME equations), which give 
a complete description of the normal and superconducting state, 
including the superconducting critical temperature. ME theory 
has been generalized to include multi-band and anisotropic cou- 
pling, magnetic and non-magnetic impurities, etc.; a review can 
be found in 1 1]. 

The biggest difficulty in the 60 's and 70 's was extracting the 
EP coupling spectral function from the available experimetal 
data, which involved some degree of approximation. In the last 
few years, it has become possible to calculate it completely ab- 
initio ||2j]. The combination of ab-initio calculations and ME 
theory has permitted to calculate the superconducting and nor- 
mal state properties of many new and old materials with consid- 
erable accuracy |3]: the biggest success is probably represented 
by the two-gap superconductor MgB2, with the record Tc of 40 
K 10, [1,0]. These methods, however, fail dramatically in more 
exotic superconductors, such as the high-r^ cuprates, where the 



key approximations (weak electronic correlations, well sepa- 
rated phonon and electronic energy scales) break down |7, sl]. 

In this paper, we analyze the possibility of applying ME 
theory to two newly-discovered pnictide superconductors, 
LaOFeAs and LaONiAs. The motivation of the application of 
the ME approach to Fe pnictides is connected with the rather 
large mass renormalisation {A^ 1 - L5 in ARPES and de Haas- 
van Alphen effects, and slightly smaller Atr ~ 0.5 in transport 
properties) as well as observed large isotope shift ape ~ 0.4 @]. 

The results for the Fe compound were already presented 
in our previous publication [10], and we review them here, 
considering also the effect of multi-band coupling. We also 
decided to include new results for the Ni compound, LaO- 
NiAs, which is the only other member of the LaOMAs fam- 
ily (M=Mn,Fe,Co,Ni,Cu,Zn) Hill showi ng supe rconductivity, 
albeit with a much lower Tc ~ 2.4-3.8 K llllfllfil. 

Experimental result s II12LI13L and calculations 

on other Ni pnictides lisl ll9|] strongly suggested that this may 
be a standard EP superconductor; in this work, we compare 
the available experimental data for LaONiAs (r^, specific heat, 
dHvA) with ME calculations, and we show that there is indeed 
a very good agreement. 

The comparison of these results with those for LaOFeAs is 
very instructive. LaOFeAs has in fact a much smaller coupling 
constant X, which is a factor 5 too low to reproduce the exper- 
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Figure 1: (color online) Crystal structure of LaOFe(Ni)As, from Ref. llOll . 



imental Tc = 26 K ll2Ql1 . even considering multiband effects. 
Similarly to the superconducting cuprates, LaOFeAs and Fe 
pnictides in general are much more "exotic" materials, where 
many-body effects may play an important role. 

The exotic features of Fe pnictides (itinerant magnetism, 
structural transitions, unusual gap symmetry, etc.) are reviewed 
in detail in other contributions to this issue, and we will not treat 
them here, although in the last part of this paper we will discuss 
how they can affect our results. 

This paper is organized as follows: In section [T]we present 
the band structure and Fermi surface of the two compounds; in 
section [2]we show the phonon dispersions and electron-phonon 
coupling calculated in Linear Response theory; in section [3] 
we present the Migdal-Eliashberg results; in the last section we 
discuss our results in light of our experimental and theoretical 
works. The technical details of the density functional theory 
(DFT) calculations are given in Appendix |A1 

1. Electronic Structure 

LaOFeAs and LaONiAs crystallize in the ZrCuSiAs struc- 
ture (space group 129); the primitive cell is tetragonal. La 
and As atoms occupy 2c Wyckoff positions, O and M atoms 
(M=Fe,Ni) occupy 2a and 2b Wyckoff positions. 

The structure, depicted in Fig. [T] consists of alternating M- 
As and La-0 layers. M and O atoms sit at the center of slightly 
distorted As and La tetrahedra; the As tetrahedra are squeezed 
in the z direction, so that there are two M-As-M angles (^i ,62), 
which are either larger or smaller than the regular tetrahedron 
value (^0 = 109.47 deg). M atoms form a square lattice; the 
M - M in-plane distance is ~ 20% larger than the M-As one. 
The relevant parameters of the structure for the two compounds 
are given in Table [T] 



In this table we report both the experimental data, from 
Ref. 1 12, 23], and the data that we obtained from a full DFT 
structural relaxation (Ref. |10] and present work), which we 
will use in the following calculations. As it was noticed by 
several authors in literature, in the Fe compound non spin- 
polarized DFT calculations tend to strongly overestimate the 
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Figure 2: Blow-up of the band structure of LaOFeAs around the Fermi level, 
decorateed with partial e (red) and t2 (green) Fe characters. From Ref. 11011 . 



As tetrahedron deformation with respect to the experiment; 
the agreement is improved if spin polarization is allowed 112 ih . 
which is normally interpreted as a sign of spin fluctuations 1221 
'23']. On the other hand, in the Ni compound, the tetrahedral 
angles given by non- spin polarized DFT calculations are very 
close to those found experimentally (see Table [TJ. 

Our band structures of LaOFeAs and LaONiAs are in very 
good agreement with literature results (l^ 21, 24, 25, 2^. The 
most important difference between the two compounds is a ~ 1 
eV shift of the Fermi level of the Ni (d^) with respect to the Fe 
(d^) compound, due to the different electron count. Measuring 
energies from the Fermi level LaOFeAs compound, O p and 
As p states form a group of 12 bands extending from ~ -6 
to -2 eV. La-/ states are found at higher energies, at ~ 2 eV. 
The dominant contribution to the states in an energy window 
extending ±2 eV around the Fermi level comes from the ten 
M-d states, which hybridize with the As p states. 

A blow-up of the band structure in this energy region, deco- 
rated with partial M d character, is shown in Figs.[2l|3l The x, y 
axes are oriented along the M - M bonds. Due to the strong 
hybridization with As p states, the d bands do not split sim- 
ply into a lower e (d^i.yi and d^^^-i) and higher t2 manifold, 
as predicted by crystal field theory. The d^i.yi orbitals, which 
lie along the M - M in-plane bonds, due to hybridization are 
split into two subsets of flat bands of located at -2 and -hi eV. 
The d^z^.i bands have the most three-dimensional character and 
sit just below Ef. The t2 bands, derived from dxy4xz^ dyz 

states, form a complicated structure centered at 0.5 eV, and 

give the largest contribution to the Density of States (DOS) at 
the Fermi level. 

The Fermi level of LaOFeAs cuts the band structure in a re- 
gion where the DOS is high (2.1 states/eV spin) and rapidly de- 
creasing; a pseudogap opens in the electronic spectrum around 
0.2 eV. The resulting Fermi surface comprises two cylindrical 
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LaOFeAs (exp) 


4.035 


8.741 


0.6512 


0.1415 


2.41 


2.85 


107.5 


113.5 


LaOFeAs (th) 


3.996 


8.636 


0.6415 


0.1413 


2.34 


2.83 


105.81 


117.1 


LaONiAs (exp) 


4.123 


8.1885 


0.6368 


0.1470 


2.35 


2.92 


103.19 


122.95 


LaONiAs (th) 


4.102 


8.2886 


0.6398 


0.1423 


2.36 


2.90 


103.99 


121.11 



Table 1: Structural data of LaOFeAs and LaONiAs from experiment (Refs. ^2^ and L12J), and DFT (Ref. (T^l and this work). Distances are in A, angles in degrees; 
for a perfect tetrahedron, ^1=^2 = 109.47 deg. 



hole pockets centered at the Y point, and a doubly-degenerate 
electron pocket centered at the M point; these sheets have a 
dominant dxz, dyz, dyz character. 

The quasi-nesting between the hole and electron pockets 
leads to a peak in the magnetic susceptibility, and hence to 
an instability of the non-magnetic solution with respect to a 
striped antiferromagnetically (AFM) ordered phase. A third 
hole pocket centered around the F point is also present; its 
character {di,^2_i or dxy) depends on the details of the calcu- 
lations, and in particular on the deformation of the As tetra- 
hedra 11221] . The plasma frequencies are strongly anisotropic 
{c4y = 2.30, ojf = 0.32 eV). 

A similar blow-up of the band structure for the Ni compound 
is shown in Fig. [3] Due to the different electron count, the xz, yz 
hole pockets are completely full, and the Fermi surface contains 
two electron and one hole sheets, with a marked 2D character. 
Besides the elliptical xz,yz pocket, the second large electron 
sheet centered at the M point has a dominant - charac- 
ter; the same bands also forms small hole pockets around the X 
point of the Brillouin zone. The corresponding bands account 
for the directional in-plane bonds of the M planes. 

The DOS is lower than in the Fe compound, A/^(0) = 
X.^^st.jeV spin, flat, and roughly particle-hole symmetric in 
an energy interval corresponding to 10 % hole and electron 
doping. The Fermi velocities are on average higher than in 
LaOFeAs, and the resulting plasma frequency are larger and 
strongly anisotropic {jufxy - 4.49 eV, ixF^ = 0.45 eV)). 

The F-centered xz, yz hole pockets are completely full, and 
this suppresses the tendency to AFM order found in the Fe com- 
pound. In fact, we do not find any AFM solution, neither in the 
LSDA nor in the GGA, in agreement with previous calcula- 
tions H. 



2. Electron-Phonon Properties 



The main results of the linear response calculations for 
LaOFeAs and LaONiAs are shown in Table |2l In both cases we 
performed the calculations at zero doping in the non-magnetic 
(NM) phase. LaONiAs is non-magnetic and superconducting 
at zero doping; on the other hand, the ground state of undoped 
LaOFeAs, both in DFT [21] and experiment, is a striped anti- 
ferromagnetic order; doping suppresses magnetism and leads 
to superconductivity. Our NM calculations are thus meant 
as a model for doped, superconducting LaOFeAs, considering 
that the eff'ect of doping in the virtual crystal approximation is 
roughly a rigid-band shift of the Fermi energy (Ef), which does 
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Figure 3: Blow-up of the band structure of LaONiAs around the Fermi level, 
decorateed with partial e (red) and t2 (green) Ni characters. 



not change the topology of the Fermi surface, but only the value 
of the DOS at E^. It is important to point out that dynamic spin 
fluctuations, which have been argued to be present also in the 
superconducting Fe samples, are not included in this calcula- 
tion, but would require going beyond the DFT level. We will 
discuss this issue in more detail in the final section of this work. 

In the two top panels of Fig. |4] we show the atom-projected 
phonon DOS (PDOS) of LaOFeAs (top) and LaONiAs (bot- 
tom). Both spectra extend up to 65 meV, and have a rather 
similar shape. The vibrations of O atoms are well separated 
in energy from those of other atomic species, lying at 61; > 40 
meV. The vibrations of La, Fe(Ni) and As occupy the same en- 
ergy range, and the eigenvectors have a strongly mixed charac- 
ter. Similarly to the electronic bands, the phonon branches (not 
shown) have very little dispersion in the z direction. Analyzing 
the evolution of the phonon eigenvectors in the Brillouin Zone 
(BZ) reveals that there is no clear separation between in and 
out-of-plane vibrations, as it often happens in layered com- 
pounds. The three major peaks in the PDOS at oj = 10,20 
and 30 meV do not show a definite in-plane or out-of-plane 
character, and cannot be easily traced back to a single vibration 
pattern. Their energy is shifted down by ~ 20 % when going 
from Fe to Ni, mostly due to EP softening. 

The EP coupling of LaONiAs (A = 0.72) is in fact much 
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2.6 


3 


BaNi2As2 


1.78 
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0.76 


4.20 
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Table 2: Electron-phonon properties of Fe and Ni superconductors calculated from Density functional perturbation theory. The results in the first two rows are from 
Ref. fid] and this work. The Tc values in the first column (T^J^) are obtained using //* = 0.12 and Allen-Dynes formula (Eq.[3}; number in parentheses correspond to 
the full numerical solution of Migdal-Eliashberg equations, given in Sect. [3] For comparison, in the last three rows we also report literature data on FeSe, LaONiP 
and BaNi2As2 from Refs. t27|,il8|,il9|]. 



larger than in LaOFeAs (A = 0.21), as shown in the two lower 
panels of Fig.lH Here we plot the two Eliashberg spectral func- 
tions a^F(<jj), together with the frequency-dependent EP cou- 
pling function A(ojy. 
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Figure 4: (color online)Frc>m top to bottom: Partial PDOS (Fe,Ni), Eliashberg 
functions and frequency-dependent electron-phonon coupling constant A(aj) 
(Eq. [TJ 12) for LaOFeAs and LaONiAs, calculated within Density Functional 
Perturbation Theory (DFPT). 
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where gv,nk,mk' are the bare DFT matrix elements, i.e. they do 
not include many-body effects. A comparison of the Eliashberg 
function with the PDOS shows that, apart from the fact that in 
both systems the high-lying O modes show very little coupling 
to electrons, there are important differences in the shape and 
size of a^F{(jS) between the Fe and the Ni compound. In fact, 
in LaOFeAs there is an almost perfect proportionality between 
the PDOS and the a^F(oj), whereas in LaONiAs the coupling 
to the two lowest peaks of the PDOS is strongly enhanced. 

A perfect proportionality between the Eliashberg function 
and the PDOS (LaOFeAs) implies that there are no patterns 
of vibration with a dramatic effect on the electronic states at the 
Fermi level. In good EP superconductors, on the other hand, the 
coupling to electrons is usually concentrated in a few selected 
phonon modes. This is best explained in terms of phonon pat- 
terns that awake dormant EP interaction between strongly di- 
rected orbitals. |28]. This is what happens in LaONiAs, where 
the electronic states derived from the Ni dx2-y2 orbitals sit at the 
Fermi level, and experience a strong coupling to the low-energy 
Ni-As modes. 

The different EP coupling of LaONiAs and LaOFeAs derives 
from the character of the electronic states at the Fermi energy. 
Therefore, it should be a rather general property of the Fe and 
Ni families of pnictide superconductors which, apart from mi- 
nor differences due to chemistry and structure, share the same 
band structure. 

In fact, as we show in Table [2l the calculated EP coupling 
constants in Ni compounds (A ^ 0.58 - 0.76) are al- 

ways 3-4 times larger than in Fe-based materials (A ^ 0.17 - 
0.21) (id Jt^Jt?]. This has important implications on the possi- 
ble pairing mechanism for superconductivity in the two classes 
of materials. We can get an estimate of Tc due to EP coupling 
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using Allen-Dynes formula 1I29I1 : 



1.2 



exp 



-1.04(1 + ^) 



^ - (1 + 0.62^)yU* 



(3) 



For yu* = 0.12, this gives Tc < 0.01 K {coin = 205 K) for 
LaOFeAs and Tc = 2.9 K for LaONiAS (coin = 96 K). 

To reproduce the experimental (Tc = 26 K) of LaOFeAs, a 
five times larger A would be needed, even for ji* = 0. On the 
other hand, the DFPT results seem to nicely explain the Tc of 
LaONiAs. 

In the next section, we will calculate in more detail the nor- 
mal and superconducting state properties of LaOFeAs and LaO- 
NiAs using the full ME theory, and compare the results with 
available experimental data. 
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3. Superconducting Properties 

a) LaFeAsOo^gFo^i 

In principle, multiband and/or anisotropic coupling could 
provide the missing factor 5 in the coupling missing to explain 
Tc similar to the multiband superconductivity in MgB2 (see, 
e.g., Refs.lS 0]); however this is very unlikely because this 
would require a very anisotropic distribution of the EP cou- 
pling llsol] . Definitely, the iron pnictides are multiband super- 
conductors with hole and electron bands which are well sep- 
arated in momentum space. To analyze this possibility we 
split the electron-phonon interaction (EPI) in Eq. [T] over elec- 
tron and hole pockets on the Fermi surface, obtaining the band- 
decomposed superconducting Eliashberg functions: 



aUco)Fij(co) = 



1 



A^/(0) 



I J,V 



(4) 



k,k',y 

S(sl)S(siMoJ-col_^.l 



where A^/(0) is the partial DOS per spin at the Fermi energy 
of the /'th sheet of the Fermi surface, is the EPI matrix 
element. These functions, shown in Fig. [51 determine the su- 
perconducting properties and thermodynamical properties like 
electronic specific heat and de Haas-van Alphen mass renor- 
malizations. In contrast to MgB2, in LaOFeAs the matrix of 
coupling constant is practically uniform: 

_ / 0.111 0.093 
~ \ 0.124 0.083 

while the characteristic logarithmic frequencies are diff'erent ( 
^intra ^ 214 K, and co\l''' = 180 K). The small diff'erence in 
the elements of the matrix A leads to a non-drastic diff'erence 
between the maximal eigenvalue and the EPI averaged over 
bands. As a result we get slightly larger value of Tc ~ 1.5 
K, for fi* = which is much lower than the observed value 
Tc = 26 K. In view of the above result, we do not pursue the 
ME study of LaOFeAs further. Other interactions, repulsive 
in the ^--wave channel but attractive in the d- or /?-wave one 
(e.g. spin-fluctuations or the direct Coulomb interaction), may 
increase Tc 1I24I1 . 
b) LaNiAsOo.gFo^i 



Figure 5: The decomposition of the EHashberg function on the intrahand and 
interband interactions; the ratio between hole and electron DOS is obtained 
from a 2"^ order fit of the band structure aroun Ep . 



LaONiAs is superconducting with Tc ~ 2.4 - 3.8 K when 
hole (Sr) or electron (F) doped. The phase diagram is roughly 
symmetric, in agreement with the flat DOS predicted by DFT 
in this doping interval. In this paragraph, we apply ME theory 
to LaNiAsOo.9Fo.i, for which we could find the most complete 
set of experimental data in literature 1 13]. 

The calculated Eliashberg spectral function a^(co)F(co) 
shown in Fig. |4l with a total A - 0.72, yields the experimen- 
tal Tc =3.8 K, with a Coulomb pseudopotential of jf = 0.12 
( slightly higher then by using the Allen-Dynes expression). 
We fix yu* = 0.12 in the following discussion. The calculated 
gap at zero temperature A(0) is 6.97 K, which gives a ratio 
2A/Tc = 3.7, higher than the BCS value. 

We now wish to investigate the temperature dependence of 
the specific heat, which yields valuable information on the size 
and nature of the EP coupling. 

In a single-band model with a strong (intermediate) EPI, in 
the normal state and in the adiabatic approximation the elec- 
tronic contribution to the specific heat is determined from the 
Eliashberg function a^(co)F(co) by the expression: 13 11I 



C'^(T) 



(2/3)n^N(0)klT 



(5) 



Jo 



l-h(6Mfir) f(col27ikBT)a^(co)F(co)co 



where N(Q) is a bare DOS per spin at the Fermi energy. The ker- 
nel f(x) is expressed in terms of the derivatives of the digamma 
function il/(x) 



f(x) = -x- 2x^5i//\x) - x^%i//''(x). 



(6) 



At low temperatures the specific heat has the well known 
asymptotic form: C^^(T 0) = (I -\- A)yoT, where A is the 
electron-phonon coupling constant, and 70 = 27i^k^N(0)/3 is 
the specific heat coefficient for noninteracting electrons. At 
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the specific heat jump we can also apply the semiempirical ex- 
pression by Carbotte 



T/T 



Figure 6: The specific heat of the LaNiASOo.gFo.i compound. The soHd(red) 
line shows results of the calculations. Black dots correspond to experimental 
results from Ref . L l3i1 



higher temperatures the specific heat diff'ers from this trivial ex- 
pression. Below Tc the diff'erence in free energies, and Fs, 
of the superconducting and normal state is given by: 



(7) 



7rA^(0)r 



2 



MZ^'iaJn) - 1) - 



^ajl(zHoJn)f+d 



where Z(ajn) is a normalization factor, (fin = An/Z((jJn) is an 
order parameter, and A„ is the gap function ( see the derivations 
in 10] and IH). 

The specific heat at temperature, T, is then calculated accord- 
ing to: 



AC,i(T) = Td\FM-Fs)/dT^ 



(8) 



The specific heat jump ACeiC^c) 8it T = Tc is determined by 
the coefficient JS = TcACqi(Tc)/2 of a second order expansion 
F^-Fs = pf, where t = (Tc - T)/Tc. 

For the comparison with experiment we have considered the 
data in Ref. llisll . The anomaly clearly visible at Tc in the 
zero-field data is suppressed by a magnetic field of 10 Tesla. 
In figure [6] the diff'erence ACp = C;,(OTesla) - Cp(lOTesla) 
is displayed as symbols. The specific heat in LaNiAsOo.Qpo.i 
was calculated using the isotropic spectral Eliashberg function 
a^(<jij)F(cL>) shown in Fig.lH The calculated specific heat at Tc 
is yN(o) = l.72ro = 6.64 mJ/molK^ with ro=3.86 mJ/molK^ 
from the band structure calculations. This value is close to the 
experimental one of 6.14 reported in Ref. lfl3l] (see Ref.Jssl]). 
The specific heat jump at Tc equals AC ^ 11.1 mJ/mol K, 
which is comparable the experimental values |13]. This gives 
AC/(y^(0)Tc) ^ 1.67 slightly larger than the BCS value of 1.43 
which corresponds to the intermediate coupling. To estimate 



Ac/r, = i.43ro(i + ^) 



l + (^)^ln(^) 

OJin 3Tc 



With ojir, = 96 K and with Tc = 3.8 K, we have AC/Tc = 
11.2 mJ/molK^, which compares well with the full numerical 
solution. The diff'erence ACei(r) = C^^(T) - C^(r) , shown in 
Figure [6] as solid line, shows a very good agreement with the 
experimental data. We would like to emphasize here that no 
fitting is involved in the theoretical calculations. The only free 
parameter which is in the Coulomb matrix element yu*, which 
was determined by the experimental Tc. 

Here we have to point out that the specific heat jumps in 
multiband systems ( or other anisotropic ones) is sufficiently 
smaller than in one-band superconductors ( see, discussions in 
1 6] and and some criteria in |35]). These results show 
that Tc and the specific heat in LaNiASOo.9Fo.i, in contrast to 
LaFeAsOo.gFo.i, can be described in the framework of the stan- 
dard single-band approach without the need of exotic mecha- 
nisms. To further support this conclusion, we observe that re- 
cent de-Haas van Alphen experiemnts on an other Ni pnictide, 
BaNi2P2, observe a band structure which is in close agreement 
with DFT calculations, with an average eff'ective mass renor- 
malization of ~ 1.8, which implies an EP coupling constant 
A = 0.8, in good agreement with our calculations [.17.1 . 

4. Discussion: 

To summarize the results of the previous sections, we have 
found that linear response calculations of the electron-phonon 
coupling yield rather diff'erent results for the two superconduct- 
ing members of the LaOMAs family of pnictide superconduc- 
tors. 

For the Fe compound the value of the total EP coupling con- 
stant A = 0.21 is much lower than in normal EP superconduc- 
tors (for example, A = 0.44 in Al, where Tc is 1.3 K), and even 
the inclusion of multiband eff'ects cannot explain the Tc=26 K 
observed in doped samples. 

For the Ni compound the coupling constant is much higher 
(A = 0.72); its normal and superconducting state properties can 
be well described by standard, single-band Migdal-Eliashberg 
theory. The values of the gap ratio 2A/ Tc = 3.7 and specific 
heat jump AC/y^(0)Tc = 1.67 are larger than what predicted 
by BCS theory (3.52 and 1.43) respectively. 

The picture that emerges from our calculations is that of 
a family of rather standard EP superconductors (Ni-based), 
opposed to a family of "exotic" superconductors (Fe-based), 
which is supported by several experimental evidences. 

The most important issue is the magnetic ground state of the 
superconducting samples. The Ni parent compounds are stan- 
dard metals, which superconduct at low T < 5 K ll^, [l3. 14. 
ISl lid llTy \j, doping (holes or electrons) does not change Tc or 



possible SDW transition has been observed at 
Ref. El 
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normal state properties dramatically. The calculated coupling 
constants /I = 0.58-0.76 |18, 19] can well explain the exper- 
imental Tc and, as we have shown in the previous section, also 
thermodynamic properties and de-Haas-van- Alphen data. 

On the other hand, in the Fe compounds superconductivity 
only appears by doping a parent compound which is an AFM 
metal. At the time when our calculations presented in Ref. fllOll 
were performed, the only available experimental data showed 
that doping suppressed the static AFM order in the supercon- 
ducting samples. |20] Therefore, we assumed that, as in the Ni 
superconductors, also in Fe superconductors the normal state is 
non-magnetic. 

In the last few months, experiments have shown that the 
ground state of the superconducting samples may also be mag- 
netic, but with fluctuating (dynamic) moments. For a more 
complete discussion of this subject, see the review by Mazin 
and Schmalian in this issue. Such an arrangement is not de- 
scribable by DFT theory, therefore it is not possible to estimate 
what its efl'ect on the EP coupling would be. 

It is also hard to compare our results with experimental data. 
There are no direct measurements of the EP coupling in liter- 
ature, although some experiments (ARPES, penetration depth, 
specifit heat) indicate some retarded electron-boson interaction, 
with a coupling constant ^0.5-1.5; however, the total cou- 
pling could be due to other bosonic excitations. 

On the other hand, a few measurements of phonon spectra are 



available in literature (M lIZl IM |39l |4Q|, |41|] . In general, there 



is a good agreement between experimental phonon frequencies 
and non-spin polarized calculations, except for the intermediate 
frequency Fe-As (and As- As) modes, which are lower in exper- 
iment than in calculations. See for example Ref. |37], where 
Inelastic Neutron Scattering data are compared to our PDOS. 
An empirical way to reconcile experiment and theory, by re- 
ducing the Fe-As force constant, was proposed in Ref. |38]. It 
was later shown that the inclusion of a static AFM order leads 
indeed to a softening the Fe-As spring constant, and improves 
the agreement of the predicted crystal structure and phonon fre- 
quencies with experiment 11421] ; however, this cannot explain the 
softening of c-polarized As modes, which form a distinct peak 
at 20 meV, at an energy 20 % lower than predicted by calcula- 
tions, which has been observed both in 1111 and 122 samples, 
and has been attributed to anomalous e - ph coupling 1391 kill. 

It was realized very early that Fe pnictide show a strong 
magneto-elastic coupling between Fe moments and As out-of- 
plane modes |25i [lOi |2l|] ; in Ref. i43l] it was proposed that this 
leads to an increased EP coupling. In principle, also many body 
efl'ects could increase the coupling constant beyond the LDA 
value. 

In conclusion, on the basis of our results we can exclude that 
standard EP coupling theory alone can cause the observed Tc in 
Fe pnictides; however, this does not mean that the phonons play 
no role in the superconducting pairing, as they might enhance 
or reduce the pairing due to other mechanisms. If this is the 
case, it is not surprising to observe a finite Fe isotope eff'ect on 

On the othere hand, LaONiAs represents a nice example of a 
single-gap, strong-coupling EP superconductor. It is important 



to stress that the diff'erence between the two compounds can be 
traced back essentially to a diff'erent filling of the same compli- 
cated, non-magnetic band structure, which derives from a non 
trivial hybridization between M and pnictogen atoms. 
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A. Computational Details: 

For the atom-projected band and DOS plots in Fig. [2ll3] 
we employed the full-potential LAPW method |44] as im- 
plemented in the Wien2k code |45]. Calculations of phonon 
spectra, EP coupling and structural relaxations were performed 
using planewaves and pseudopotentials with QUANTUM- 
ESPRESso 1 46]. We employed ultrasoft Vanderbilt pseudopoten- 
tials 1 47], with a cut-off' of 40 Ryd for the wave-functions, and 
320 Ryd for the charge densities. The k- space integration for 
the electrons was approximated by a summation over a 8 8 4 
uniform grid in reciprocal space, with a Gaussian smearing of 
0.02 Ryd for self-consistent cycles and relaxations; a much finer 
(16 16 8) grid was used for evaluating DOS and EP linewidths. 
Dynamical matrices and EP linewidths were calculated on a 
uniform 442 grid in q- space; phonon dispersions and DOS were 
then obtained by Fourier interpolation of the dynamical matri- 
ces, and the Eliashberg function by summing over individual 
linewidths and phonons. To check the eff'ect of nesting on the 
EP linewidhts, we also calculated selected q points on a 882 
grid. 

Whenever possible, we cross-checked the results given by 
the two codes and found them to be in close agreement; for 
consistency, we used the same GGA-PBE exchange-correlation 
potential in both cases ll48ll . 
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